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This study investigates the dynamic relationship between coal consumption, economic growth, trade 
openness and CO, emissions in case of India. In doing so, Narayan and Popp, Journal of Applied Statistics 
2010; 37:1425-1438, structural break unit test is applied to test the order of integration of the variables. 
Long run relationship between the variables is tested by applying the ARDL bounds testing approach to 
cointegration developed by Pesaran et al. Journal of Applied Econometrics 2001; 16:289-326. 

The results confirm the existence of cointegration for long run between coal consumption, economic 
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Energy mental Kuznets curve (EKC) in long run as well as in short run. Coal consumption as well as trade 
eae openness contributes to CO3 emissions. The causality analysis reports the feedback hypothesis between 
EKC economic growth and CO, emissions and same inference is drawn between coal consumption and CO2 
emissions. Moreover, trade openness Granger causes economic growth, coal consumption and CO3 

emissions. 
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1. Introduction 


The Kyoto Protocol requires that industrialized countries 
reduce their collective emissions of greenhouse gasses by 5.2% 
of 1990 levels by the period 2008-2012. The country-specific 
targets in the Kyoto Protocol may be difficult for some nations to 
achieve. Developing countries, including India, have absolved of 
any responsibility towards reducing emissions in the first com- 
mitment period, that is, 2008-2012, of the Kyoto Protocol. 
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However, India estimates a 24% cut by 2020 in its carbon intensity 
i.e., the amount of carbon dioxide emitted for each unit of GDP, 
compared with 2005 levels, and by 2030, it estimates it could 
achieve a reduction in its carbon emissions by 37% from 2005 
levels, according to provisional government figures (The Indian 
Express [71]). 

Theoretically, the environmental Kuznets curve (EKC) hypothesis 
postulates the existence of an inverted-U shape relationship 
between real GDP per capita and measures of environmental 
degradation such as SO2 and/or CO% emissions. However, the 
empirical evidences, either using time series and/or pooled data of 
a group of countries, on the EKC hypothesis vary from country to 
country, instead. Further, the results are not uniform across pollu- 
tants. This has created tow problems being faced by environmental 
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policy makers: to ensure that useful knowledge informs policy 
(without being misused and/or distorted) and to understand how 
to respond to this knowledge (Boehmer-Christiansen [1]). However, 
in the present context we have limited ourselves to provide 
evidence of the EKC hypothesis for India. Our contribution lies, 
particularly in two directions. First, we employed the ARDL method 
which is amenable for short time series data as in this paper and 
second, we provide empirical evidence of the EKC by incorporating 
coal consumption and trade as additional determinant of CO2 
emissions in case of India. 

The rational for selecting India for our analysis is that it has 
implemented a variety of programs and policy initiatives since 
the introduction of the National Forest Policy in 1988. India has 
become one of the fast growing countries next to China. However, 
in the same time consumption of coal has increased rapidly which 
was 35.55 mote in 1965 has reached the height of 249.86 mote in 
2009. The close relationship between coal consumption and 
economic growth of India is evident from following Fig. 1. Further 
if we see the relationship between percentage growth rates of 
GDP and coal consumption, we find that recently percentage 
growth rate in the coal consumption is more than twice of the 
percentage GDP growth rate. This situation is evident in Fig. 2. 

With this background we set objective to test the environ- 
mental Kuznets curve hypothesis (EKC) in augmented equation 
framework. According to the EKC hypothesis environmental 
degradation increases at initial level of economic growth and 
then starts to decrease at a higher level of economic growth. 
Hence, the relationship between measure of environmental 
degradation, (in our case it is CO emissions) and measures of 
economic growth (in our case it is measured by real GDP per 
capita) is the inverted-U shaped curve. The objective of present 
study is to investigate the EKC for the Indian economy over the 


Coal 


a GDP per capita 


Fig. 1. Coal consumption and GDP per capita. 
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Fig. 2. Percentage growth rate in coal consumption and GDP. 


period of 1966-2011. In addition, we also include coal consump- 
tion and trade openness. 

The rest of the paper consists of: Section 2 reviews a selected 
literature encompassing the EKC and the variables listed above, 
Section 3 has the theoretical and the econometric model includ- 
ing the ARDL estimation strategy. The empirical results are 
reported in Section 4 followed by the conclusion and policy 
implications. 


2. Literature review 


There has been increasing attention on the impact of economic 
growth on environment degradation since the last few decades of 
the previous century. The pioneering attempt in this area was 
made in the early 1990s by Grossman and Krueger [2,3] which 
investigated the environmental impacts of the North American 
Free Trade Agreement. Their study postulated, estimated and 
ascertained an inverted-U shaped relationship between measures 
of several pollutants and real GDP per capita which was, con- 
temporaneously, confirmed by Shafik and Bandyopadhyay [4] and 
Panayotou [5]. Theoretically, the EKC does not only depend on 
levels of per capita GDP, but also on a series of factors causing 
changes in economic growth that can affect environment. In 
general, economists analyze mechanisms behind the EKC by 
examining scale effect, structural effect and technique effect 
(Song et al. [6]). However, few scholars (for example, Panayotou 
[5]) believe that EKC is caused by up-gradation from the adjust- 
ment of economic structure. According to the structural effect 
hypothesis (Stern [7]) economic development passes through the 
various states starting form preliminary, rapid-development and 
high-grade, industrial structure first upgrades from agriculture to 
a high-pollution industry, finally turns to information concen- 
trated industry, which leads to the improvements in environ- 
mental quality. However, as Stokey [8] pointed out, due to the 
technique effect economic growth can break through one thresh- 
old point after arriving at a certain stage of economic develop- 
ment. Hence, at a low income level, only the high pollution 
technique can be used, but once leaping over the threshold point 
of economic development, cleaner technologies can be adopted 
which lowers the degradation in the environmental quality. 
Further, some scholars attribute the demand factors to the cause 
of EKC (for example, Lopez [9]), which asserts that demand for a 
clean environment will be increased over the real income per 
capita. Giving the importance of scale effect, Andreoni and 
Levinson [10] suggested that in the static model of single depart- 
ment, the EKC can be derived technically, only if pollution control 
is increasing in scale. 

Further, Suri and Chapman [11] bring the contribution of 
industrial products of the imports and exports to industrial 
products of national production into the analytical framework 
of the EKC. That is to say, the low emissions correspond to the 
growth of industrial products of the imports, while the high 
emissions correspond to that of the exports. This implies that 
there is a strong relationship between trade and environmental 
quality, so the evolution of environmental quality can be pre- 
dicted effectively. The empirical studies of EKC started by Gross- 
man and Krueger [2] and followed by Lucas et al. [12], Wyckoff 
and Roop [13], Suri and Chapman [11], Heil and Selden [14], Friedl 
and Getzner [15], Stern [7], Nohman and Antrobus [16], Dinda 
and Coondoo [17] and Coondoo and Dinda [18] but presented 
mixed empirical evidence on the validity of EKC. Song et al. [6], 
Dhakal [19], Jalil and Mahmud [20] and, Zhang and Cheng [21] 
supported the existence of EKC in China. The findings of Fodha 
and Zaghdoud [22] revealed the existence of EKC between the SO, 
emissions and economic growth but not for the CO, emissions in 
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Tunisia. In contrast, Akbostancl et al. [23] did not support the 
existence of EKC in Turkey and same conclusion was drawn by 
Saboori et al. [64] for Indonesia. They argued that CO, emissions 
are automatically reduced due to the rapid pace of economic 
growth. 

Recent literature documented alliance of economic growth 
with energy consumption and environmental pollution to inves- 
tigate the validity of EKC. The relationship between economic 
growth, energy consumption and COz emissions have also been 
researched extensively both in the country case and panel studies. 
For example, Ang [24] found stable long run relationship between 
economic growth, energy consumption and CO, emissions for 
French economy while Ang [25] also got similar results for 
Malaysia. Ang [24] showed that causality is running from eco- 
nomic growth to energy consumption and CO, emissions in the 
long run but energy consumption Granger-causes economic 
growth in the short run. In case of Malaysia, Ang [25] reported 
that output increases CO2 emissions and energy consumption. 
Ghosh [26] documented that no long run causality between 
economic growth and CO, emissions and in the short run, 
bidirectional causality exists in India. On contrary, Tiwari [66] 
noted unidirectional causality running from CO, emissions to 
economic growth and CO2 emissions and energy consumption 
and CO, emissions Granger cause each other. In case of Indonesia, 
Hwang and Yoo [69] noted bidirectional causality between energy 
consumption and CO, emissions in long run but in short run, 
energy consumption and CO, emissions are cause of economic 
growth. Latter on, Uddin et al. [70] investigated the relationship 
between energy consumption, economic growth, trade openness 
and CO; emissions in case of Srilanka. Their results found that 
economic growth Granger causes energy consumption and CO, 
emissions. 

For the panel studies, Apergis and Payne [27] investigated this 
relationship for six Central American economies using panel 
VECM. It is evident that energy consumption is positively linked 
with CO, emissions and EKC hypothesis has been confirmed. Lean 
and Smyth [28] and Apergis and Payne [29] reached the same 
conclusion for the case of ASEAN countries and Commonwealth of 
Independent States, respectively. Narayan and Narayan’s [30] 
empirical evidence also validates the EKC hypothesis for 43 low 
income countries. In addition, Lean and Smyth [28] noted a long 
run causality running from energy consumption and CO emis- 
sions to economic growth but in short span of time, energy 
consumption Granger causes CO, emissions. On the other hand, 
Apergis and Payne [29] found that energy consumption and 
economic growth Granger causes CO, emissions while bidirec- 
tional causality is found between energy consumption and 
economic growth; and between energy consumption and CO, 
emissions. Chen [31] explored this issue to Chinese provinces and 
documented that industrial sector’s development is linked with 
increase of CO, emissions due to energy consumption. Similarly, 
Shiyi [32] investigated the relationship between industrial sec- 
tor’s development and CO, emissions using Chinese provincial 
data and concluded that industrial development increases CO2 
emissions. Pao and Tsai [33], Ozturk and Acaravci [34] and, 
Acaravci and Ozturk [35] also validated the existence of the EKC 
in case of BRIC, Turkey and Demark and Italy, respectively. In case 
of Pakistan, Nasir and Rehman [36] and, Shahbaz et al. [37] 
reported that the EKC exists while a rise in economic growth 
and energy consumption is linked with an increase in energy 
pollutants. Iwata [38] investigated the empirical existence of the 
EKC in 28 countries adding Pakistan. They corroborated the 
existence of the EKC by nuclear energy consumption taking into 
account. Saboori et al. [65] unveiled that long run as well as short 
run EKC exists in case of Malasia. Unidirectional causality from 
economic growth to CO z emissions further malidates the 


existence of KEC for Malaysian economy. Yeh [67] used the 
quintiles regressions using the data of developed and developed 
countries and reported the existence of EKC by controlling other 
macroeconomic variables. In case of Romania, Shahbaz et al. [68] 
confirmed long run relationship between economic growth, 
energy consumption and energy pollutants. Their empirical evi- 
dence also found that environmental Kuznets curve (EKC) exist 
both in long-and-short runs. Moreover, energy consumption is 
major contributor to energy pollutants. Democratic regime shows 
her significant contribution to decline CO, emissions through 
effective implementation of economic policies and financial 
development improves environment i.e., reduces CO, emissions 
by redirecting the resources to environment friendly projects. 

The relationship between trade openness and environmental 
degradation has also been investigated empirically. Grossman and 
Krueger [2]) argued that environmental effect of international 
trade depends on the policies of an economy. There are two 
schools of thought about the impact of international trade on CO2 
emissions. One argued that trade openness provides an offer to 
each country to have access to international market which 
enhances the market share among countries. This leads the 
competition among the countries and increases the efficiency of 
using scarce resources and encourages importing cleaner tech- 
nologies to decline COz emissions (e.g., Runge [39] and Helpman 
[40]). Other group probed that natural resources are depleted due 
to international trade. This depletion of natural resources raises 
CO» emissions and causes environment quality worsened (e.g., 
Schmalensee et al. [41]; Copeland and Taylor [42]; Chaudhuri and 
Pfaff [43]). In country case studies, Machado [44] indicated 
positive link between foreign trade and CO emissions in Brazil. 
Mongelli et al. [45] concluded that pollution haven hypothesis is 
existed in Italy. Halicioglu [46] added trade openness to explore 
the relationship between economic growth, CO emissions and 
energy consumption for Turkey. The result showed that trade 
openness is one of main contributor to economic growth while 
income raises the levels of CO emissions. Finally, Jayanthaku- 
maran et al. [47] probed the relationship between economic 
growth, energy consumption, trade openness and CO2 emissions 
in case of India and China. Their empirical evidence confirmed the 
validation of EKC in both countries. 


3. Modeling, empirical strategy and data collection 


The dynamic relationship between income, coal consumption 
and CO; emissions is investigated by including trade openness as 
potential determinant to CO, emissions. For this purpose, we 
followed the empirical model applied by Ang [25,26] for France 
and Malaysia, Soytas et al. [48] for United States, Jalil and 
Mahmud [20] for China, Halicioglu [46] for Turkey, and Shahbaz 
et al. [37] for Pakistan. We use coal consumption as an indicator 
of energy rather than energy consumption. More than 50% energy 
demand in India is fulfilled from coal consumption. Following 
Shahbaz et al. [37], we converted all the series into natural 
logarithms to obtain efficient and consistent results. The relation 
is specified as follows: 


CO2 =f (Ye Y?, Ce TR) (1) 
In the log-linear form the specification is written as: 


InCO.; = h1 + 7InY;+ @3lnY? + b4lnC; + dsINTR; + he (2) 


where, CO2 denotes carbon emissions per capita (in kt); C is coal 
consumption per capita; Y and (Y?) mentions real GDP per capita 
and its square, respectively. TR is trade openness [(exports+im- 
ports)/GDP)] per capita u is a random error term. Following the 
EKC hypothesis we expect that p, > 0 and p; <0. The use of coal 
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is detrimental for environment and it is expected that $, > 0. The 
impact of trade openness on CO, emissions can be positive or 
negative. This implies that ¢, <0 if government enforces envir- 
onmental laws, and makes possibility to import of environment 
friendly capital and technology to be used in production process. 
But Grossman and Krueger [2,3] argued that if emissions might be 
generated due to relocation of polluting industries from devel- 
oped economies, a practice known in the literature as the ‘safe- 
haven hypotheses’ the ¢; > 0. 

We applied a recent unit root test developed by Narayan 
and Popp [49] to attain efficient and consistent empirical evi- 
dence about the integration order of the variables. This test 
incorporates about two structural breaks stemming in the series 
which is based on two regression equations of the form as 
follows: 


yi = py, +o tp t+ O1D(Tg),,+92D(Tp) 9. 


k 
+ô1DU; t1 +82DU, t1 + X BAY +e (3) 
j=l 


YE? = PY 1 +*+ B*t+11D(Tp), + K2D(Tp) 2, 
+0;DU; t1 +65DU5 t1 + YTDT) 1 +73DT 2-1 


k 
+ So BAY j+ee (4) 
j=l 


where e, is a white noise process, such that es~ NIID(0,07), 
D(Ts) + =1(t=T3;+1), i=1,2, DU; =1(t>Tgi) DTi,= 
1(t > T3;)(t—DT3;), Tgi i=1, 2, denote the true break dates. Eqs. 
(3) and (4) are IO-type test regression for M1 and M2, respec- 
tively. The unit root null hypothesis of p=1 is tested against the 
alternative hypothesis of p < 1 by using the t-statistics of 6. Since 
it is assumed that true break dates are unknown, Tg; in Eqs. 
(3) and (4) has to be substituted by their estimates Tai, i=1, 2, in 
order to conduct the unit root test. We used a sequential 
procedure following Narayan and Popp [49] to attain Tai, i=1, 2. 

We apply the autoregressive distributed lag model or the ARDL 
bounds testing approach to cointegration. The ARDL bounds 
testing approach to cointegration is preferred due to its certain 
advantages over traditional cointegration approaches. For exam- 
ple, the ARDL bounds testing is flexible regarding the integrating 
order of the variables whether variables are found to be stationary 
at I(1) or I(0) or 1(1)/I(0)'. The Monte Carlo investigation shows 
that this approach is superior and provides consistent results for 
small sample (Pesaran and Shin [50]). Moreover, a dynamic 
unrestricted error correction model (UECM) can be derived from 
the ARDL bounds testing through a simple linear transformation. 
The UECM integrates the short run dynamics with the long run 
equilibrium without losing any information for long run. The 
empirical formulation of the ARDL bounds testing approach to 
cointegration is given below: 


AlnCO? = 041 +arT + aco INCO? ; 
+ dy InY+_4 + y2 Iny?., + acInC;_4 + Orrin TR: 


p q r 
+ X ajAlncO?_,+ XO ajAlnY,_; + XO %AlnY? , 
j=0 


i=1 k=0 
Ss t 
+ 5 oAlnC,_; + 5 AmAln TRe-m + he (5) 
1=0 m=0 


AlnY; = 01 +4rT+&co INCO? 4 + ayINY;-1 


1 The ARDL approach to cointegration is applicable if variables are integrated 
at I(1) or I(O) or I(1)/I(0). 


p 
+ OyalnY?_, +occlnCr_1 +oereln TRea + XO BAlnY;_j 
= 


q r S 
+ XO BjAlncoz_, + XO B,AlnY¢_, + XO BAlnC,_; 
j=0 k=0 l=0 


t 
+ > BmAln TR, + My (6) 


m=0 
AlnY? = %1 + &rT + &¢q2INCO?_ + ayInY;_1 


p 
+aylnY? ; +aclnCe1 +ærrln TR + D B,AlnY2_, 
=l 


q r s 
+ XO BjAlInCOz_,+ S> BeAINY;_~ + X BANC 
j=0 l=0 


k=0 


t 
+ 5 Bm AlN TRe-m + My (7) 


m=0 
AlnC; = %1 +7T + %¢g2InCO?_, +ayInY +1 


p 
+ Ay2 Iny?_, + acInC;_4 + Orpln TR- + 5 Q;AlnC;i 
i=1 


q r š 
+ $O 9Alncoz_, + XO 9, AlnY,_~ + X 9AlnY? 
j=0 1=0 


k=0 


t 
+ JO SmAln TRem+ bp (8) 


m=0 
Aln TR; = 01 + &rT + O¢¢q2InCO?_, +ayInY;_1 4+ %y2InY?_, 


p 
+aclnCe + oyeln TR + 5 p;^ln TRei 
i=1 


q r s 
+ X p,Alnco?_, + XO p,AlnY,_~+ X p,AlnY?_, 
j=0 k=0 i=0 


t 
+ JO PmAINCrm+ Me (9) 
m=0 


here, to test for the existence of cointegration calculated F- 
statistic is used to compare with critical bounds generated by 
Pesaran et al. [51]. Two critical bounds namely upper critical 
bound (UCB) and lower critical bound (LCB) have been developed 
by Pesaran et al. [51]. In our specification of Eqs. (5)-(9), the null 
hypothesis which is tested i.e., Ho : & gz = Xy = y2 = Ac = Arr = 0 
against alternate hypothesis Ho : &o2 # Hy A y2 # Ac # Arr # 0 to 
test for cointegration. The F-test is non-standard and we may use 
the LCB and UCB developed by Pesaran et al. [51]? to test for 
cointegration. If computed F-statistic is more than UCB, there is 
cointegration between the variables. If computed F-statistic does 
not exceed the LCB, there is no cointegration relationship and if 
computed F-statistic falls between lower and upper critical 
bounds then decision regarding cointegration between the vari- 
ables remains uncertain? . However, the critical bounds generated 
by Pesaran et al. [51] may not be appropriate for small sample like 
ours case with 45 observations hence, we used critical bounds 
developed by Narayan [52]. Further, we have conducted the 
stability tests to examine the stability of ARDL bounds testing 
estimates by utilizing CUSUM and CUSUMSQ test of (Browns et al. 
[53]). 


2 Pesaran et al. [51] have computed two asymptotic critical values—one when 
the variables are assumed to be I(0) and the other when the variables are assumed 
to be I(1). 

3 In such case, error correction method is appropriate method to investigate 
the cointegration (Banerjee et al. [54]). This indicates that error correction term 
will be a useful way of establishing cointegration between the variables. 
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After confirming the evidence of cointegration and cointegra- 
tion is stable we moved ahead for testing the causal relationship 
between the test variables in the framework of vector error 
correction method (VECM). The VECM in five variables case can 
be written as follows: 


l m 
Aln co? = 001+ 5 4, Aln co? + X &z Aln Yı-j 


i=] j=l 


n o 
+ > 033 Aln Ya + 5 ota4AlnC;_, 


k=1 r=1 


p 
+ $5 as5Aln TRLi-s +19, ECT ra + Hi (10) 


s=1 


l m 
AlnY = Bo, + X` By, Aln Y;_j+ XO Bp2Aln COF 


i=1 j=l 
n o P 
+ 52 BsgAln Y2 4+ XO paa AnCir+ XO BssAlnTR,_1 
k=1 r= s=1 
+ NQECT 1 + Hi (11) 


l m 
Aln Y? = $o + ` n Aln YẸ + D $22Aln COP, 


i=1 j=1 
n o P 
+ 5 h33AIn Yek + 5 Pa, AlnCi1 + > Ps5AlNTR:_1 
k=1 r=1 s=1 
+3ECT;_1 + Ls; (12) 


I m n 
AlnC: = $01 + X dy AlnCi1+ X 22AlnCOz_; + XO G33AInY; 
i=1 j=l k=1 


o p 
+5 Pa4AlnY? 1 + XO PssAINTRe 1 +N gECT ra + Mai 


r=1 $= 
a3) 
l m 
AInTR; = 691 + 5 61, AlnTR,_;+ 5 ô22AlnCO2 ; 
i=l jaza 
n o P 
+ 5 33 AlNY +k + 5 ô44AlnY? + 5 dssAlnC;_4 
k=1 r=1 s=1 
+gECTi 4 + Maj (14) 


where u; are white noise residual terms. The estimates of ECT;_1 
also shows the speed of convergence from short run towards long 
run equilibrium path in all models depending upon the sign of the 
coefficient of ECT,.; term. The usefulness of the VECM lies in 
distinguishing between short-and-long runs causal relationships 
in one hand and detecting causality in long-run, short-run and 
joint i.e., short-and-long runs on the other hand. The significance 
of ECT;_; can be tested using t-statistic which confirms the long- 
run causal relationship whereas the short run causality is 
detected by the joint x? statistical significance of the estimates 
of first difference lagged independent variables. In the final step 
we used Wald-test to test the joint significance of estimates of the 
lagged terms of independent variables and error correction term 
which confirms the existence of short-and-long run causality 
relations (Shahbaz et al. [55]) and known as measure of strong 
Granger-causality (Oh and Lee [56]). 

The data on these variables has been collected from World 
Development Indicators (CD-ROM, 2011). Further, World Devel- 
opment Indicates has been used to collect data on real GDP per 
capita, coal consumption per capita, real trade (real imports per 
capita+ real exports per capita) per capita and CO, emissions per 
capita. The study covers the period of 1966-2009. 


Table 1 
Ng-Perron unit root test. 


Level 
Variables MZa MZt MSB MPT 

InCOz; —9.8577 —2.2095 0.2241 9.2910 
InY; — 8.9336 — 1.8632 0.2085 11.0896 
Iny2 1.2003 0.6869 0.5722 85.0399 
Inc; — 4.1594 — 1.1484 0.2761 19.0857 
InTR: — 4.8699 — 1.4368 0.2950 18.0242 
1st difference 

InCO2; —24.7353* —3.5163 0.1421 3.6862 
InY; —19.2482** —3.1020 0.1611 4.7358 
Iny? —23.9443* — 3.4600 0.1445 3.8058 
Inc; —21.5425** —3.2815 0.1523 4.2324 
InTR; —22,4313* —3.3472 0.1492 4.0729 


* Indicate the significance at the 1% level. 
** Indicate the significance at the 5% level. 


4. Empirical results and discussions 


We have applied the ARDL bounds testing approach to coin- 
tegration to find long run relationship between economic growth, 
coal consumption, trade openness and CO, emissions in case of 
India. Although bounds testing approach is flexible regarding the 
integrating order of the variables. We can apply it, if the variables 
are integrated at I(0), I(1) or I(0)/I(1). We cannot compute the 
ARDL F-statistics for long run relationship if any variable is found 
to be stationary at I(2). To ensure that none of the variables is 
stationary at I(2) or beyond that order of integration, we applied 
Ng-Perron [57] unit root test. The results are reported in Table 1. 
This empirical evidence indicates that all the series have unit root 
at their level form but found to be stationary at 1st difference. 
This shows that all the variables have I(1) order of integration. 
The empirical evidence provided by Ng-Perron may be biased 
because it does not have information about structural break 
stemming in the series. 

To solve this issue, we applied unit root test with two 
structural breaks developed by Narayan and Popp [49]. The 
inappropriate chose of structural breaks may reject null hypoth- 
esis and provide inconclusive results. The technique used in N-P 
unit root test allows in determining two structural breaks in the 
level and trend of the series endogenously. The unit root test 
developed by Narayan and Popp [49] indicates two structural 
breaks efficiently as compared to other structural break unit root 
tests. Table 2 reported the results by Narayan and Popp [49] unit 
root test. Our empirical exercise reveals that all the variables are 
found to be non-stationary (except InTR;, which is in model 
2 turns to be stationary) at level but integrated at 1st difference. 
This confirms that all the variables are integrated at I(1) and 
robust‘. 

Our empirical evidence confirms that all the variables are 
stationary at 1st difference form that intends to apply the ARDL 
bounds approach in examining the long run relationship between 
the series. It is prerequisite to choose appropriate lag length of the 
variables to capture the dynamic relationship between the vari- 
ables. The statistics of F-test are very much elusive with the 
selection of lag order of the variables. Inappropriate lag length 
selection provides biased results and in resulting, decision about 
cointegration is worthless (Ltitkepohl [58]). We have used AIC and 
SBC criterion to choose lag order but our decision is based 
minimum value of AIC. Litkepohl [58] suggested that lag order 
selected by AIC provides efficient and consistent results. The 


4 The results of N-P at 1st difference are not provided here due lack of space 
and are available from authors upon request. 
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Table 2 
Narayan and Popp [49] Structural Break Unit Root Test. 


Table 5 
Long-and-short run analysis. 


Variable Model 1(M1) Model 2(M2) 

T-statistic TB1 TB2 T-statistic TB1 TB2 
InCO3,; —4.0378 1976 2000 — 4.3381 1976 2000 
InY; 0.7720 1974 1978 0.4483 1974 1978 
Iny? 1.0668 1974 1978 0.7999 1974 1978 
Inc; —1.1955 1973 1989 —1.6512 1977 1989 
InTR; —3.0225 1974 1981 — 5.3848” 1984 1994 


Note: Model-1 assumes two breaks in level and Model-2 assumes two breaks in 
level as well as slope. 


** Denotes significance at 5% level. 


Table 3 
Results of ARDL cointegration test. 


Variable InCO, Iny; Iny2 Inc; InTR; 
F-statistics 8.9129% 7.1597" = 7.1578"** 7.4063™™ 5.3858 
Critical values* 1% level 5% level 10% level* 

Lower bounds 10.150 7.135 5.950 

Upper bounds 11.130 7.980 6.680 

Diagnostic tests 

R2 0.8325 0.9994 0.9995 0.6249 0.8018 
Adj—R? 0.6802 0.9990 0.9991 0.3796 0.6037 
Durbin-Watson 1.6951 2.2267 2.2155 2.0707 2.3233 


* Critical values bounds are from Narayan [52] with unrestricted intercept and 
unrestricted trend. 

* Show the significance at 1% level. 

** Show the significance at 5% level. 

*** Show the significance at 10% level. 


Table 4 
Results of Johansen cointegration test. 


Hypothesis Trace statistic Maximum eigen value 
R=0 105.8022" 46.7978" 

R<1 59.0044" 32.3417" 

R<2 26.6626 17.6688 

R<3 8.9938 8.6642 

R<4 0.3295 0.3295 


* Show significant at 1% level. 
** Show significant at 5% level. 


reason is that AIC has superior properties and high explaining 
power as compared to SBC. Our empirical evidence shows that we 
cannot take lag more than 2 because data sample consists of 45 
years Table 3. We have used yearly data over the period of 1966- 
2009 of Indian economy”. 

After choosing lag order, next step is to compare our computed 
F-statistics with lower and upper critical bounds generated by 
Narayan [52]. The results show that our calculated F-statistics 
exceed upper critical bounds at 5% and 10% levels, respectively, 
once, COz emissions per capita, real GDP per capita (squared of 
real GDP per capita) and coal consumption per capita treated as 
respond variables. The results are reported in Table 4. 

This implies that three cointegrating vectors are found recom- 
mending the long run relationship between coal consumption, 
economic growth, trade openness and CO, emissions in case of 
India over the study period of 1966-2011. 

The robustness of long run results established by the ARDL 
bounds testing approach to cointegration is tested by applying 


3 Results of lag length selection are available from authors upon request. 


Dependent variable =!nCO,, 


Long run results 


Variable Coefficient Std. Error T-Statistic 
Constant —50.2450 4.7539 —10.5692* 
InY; 10.0674 0.8740 11.5182* 
Iny? — 0.4942 0.0410 —12.0467* 
nC: 0.8436 0.0999 8.43771* 
InTR 0.0862 0.0415 2.07734** 
Short run results 

Constant 0.0282 0.0055 5.0598* 
ny; 4.9768 1.6149 3.0816* 
ny? — 0.2392 0.0807 —2.9622* 
InC; 0.4481 0.0948 4.7244* 
nT Re —0.1090 0.0453 —2.4035** 
ECM:-1 —0.2994 0.0948 —3,1574* 
Diagnostic tests 

Test F-statistic Prob. value 

X2 NORMAL 0.8323 0.6595 

X SERIAL 1.9906 0.1548 

y2 ARCH 0.0150 0.9030 

y2WHITE 0.7055 0.7110 

y? REMSAY 2.3963 0.1321 


* Denote the significant at 1% level. 
** Denote the significant at 5% level. 


Johansen multivariate cointegration approach. The results 
reported in Table 4 confirm that there are two cointegrating 
vectors at 1% level of significance indicated both by Trace statistic 
and Maximum Eigen values. This corroborates the long run 
relationship between coal consumption, economic growth, trade 
openness and CO, emissions in case of India over the study period 
of 1966-2011. 

After investigating the robustness of long run relationship 
between these variables, the long run marginal impacts of coal 
consumption, economic growth, trade openness on CO, emissions 
in case of India is investigated by applying OLS approach. The 
results reported in Table 5 postulated that coal consumption is 
major contributor to energy pollutants and it is statistically 
significant at 1% level of significance. A 1% rise in CO2 emissions 
is associated with 0.8436% consumption of coal. These findings 
support to the view by Wolde-Rufael [59] that coal consumption 
highly contributes to energy pollutants as compared to other 
indictors of energy. Jayanthakumaran et al. [47] also nominated 
that India meets 55% of energy demand through coal consump- 
tion and produces 5% of CO, emissions of the world®. Further- 
more, Jayanthakumaran et al. [47] reported that in total energy 
emissions 86% is contributed by coal consumption in case of India 
and same inference drawn by Alam et al. [60] . 

The impact of trade openness on energy emissions is positive 
and significant at 5% level. All else is same, a 1% expansion in 
trade openness leads to upsurge in energy emissions by 0.0862%. 
This finding is consistent with Jayanthakumaran et al. [47] also 
reported positive effect of international trade on CO2 emissions 
but it is statistically insignificant. In regional studies, Halicioglu 
[46] reported that international trade increases energy emissions 
by raising economic growth and same inference is drawn by Nasir 
and Rehman [36] in case of Turkey and Pakistan, respectively. On 
contrary, Shahbaz et al. [37] noted that trade openness improves 
the environmental quality by enhancing the capacity of country 


6 Other sources of energy such as crude oil, natural gas and petroleum 
products are used to meet the rest 45% energy demand. 
7 Alam et al. [60] used total primary consumption as an indicator of energy. 
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Table 6 
VECM Granger causality analysis. 


Dependent variables Direction of causality 


Short run Long run Joint (Long-and-short runs) causality 
InCO2,; InY; Iny? Inc; InTR ECTy_4 InCO2,,ECT;_, InY;,ECT,_1 InY?,ECT;_1 InC;,ECT,_1 INTR, ECT: 
nCOx = 3.6272" 3.3713** 12.8363" 5.1278** —0.4550* = - 7.5935* 6.9193* 9.3793* 7.8799* 
0.0384] [0.0473] [0.0001] [0.0119] [—3.9339] [0.0006] [0.0040] [0.0001] [0.0005] 
<mmEquation id="eqn87” | 4.6517** 8.1920** 1.5550 6.1098* —0.1632** 3.4984** = 5.8590* 1.2920 4.5041** 
> [0.0174] [0.0000] [0.2227] [0.0060] [—1.8427] [0.0274] [0.0019] [0.2952] [0.0101] 
nY? 4.2919* 8.4460* - 0.9678 5.8392* —0.1535*** 3.1922** 6.0750* - 1.1602 4.3125" 
[0.432] [0.0001] 0.3878] [0.0072] [—1.7716] [0.0377] [0.0045] [0.3412] [0.0121] 
nC 12.5233" 0.5701 0.5177 - 2.7016* —0.7704* — 10.0023* 6.5615* 6.6326* z 7.2152* 
[0.0001] [0.5715] [0.6013 [0.0834] [—4.2526] [0.0001] [0.0015] [0.0012] [0.0009] 
nT re 2.4178" 4.6856** 4.5034** 1.0847 - 2 z = - = = 
[0.1078] [0.0169] [0.0195] [0.3509] 
* Show significant at 1% level. 
** Show significant at 5% level. 
*** Show significant at 10%,.level in parenthesis we reported p-values. 
to implement advanced technology to increase domestic ——— CUSUM ——— 5% Significance 
production. ee Le 
The impact of linear and nonlinear terms of real GDP per capita ‘ We estes a 
is positive and negative. The coefficients of linear and nonlinear 4 Leseneeseeneese 
terms are 10.0674 and — 0.4942 and it is statistically significant at 0 
1% level of significance, respectively. The significance of both 44 
linear and non-linear terms of real GDP per capita provides the ce ita ee 
empirical evidence of inverted-U shaped relationship between MNOS en a 
economic growth and CO% emissions “so called environmental aG I 1990 1905 000 005 mn 


Kuznets curve” (EKC). The empirical exercise revealed that a 
10.0674% energy emissions are linked with 1% increase in GDP 
per capita and inverse effect of squared term of real GDP per 
capita shows the delinking point of CO, emissions i.e., — 0.4942 
and real GDP per capita, once an economy achieves higher level of 
income. This empirical evidence provides the support for EKC 
revealing that CO emissions increase at the initial stage of 
economic growth and decline after a threshold point i.e., Ind. 
Rupee 28,131 (US$531). These results are with existing literature 
such as He [61], Song et al. [6], Jalil and Mahmud [20], Fodha and 
Zaghdoud [22] and Jayanthakumaran et al. [47] for China and 
India, Halicioglu [46] for Turkey, Lean and Smyth [28] for ASEAN 
countries, Nasir and Rehman [36] and Shahbaz et al. [37] for 
Pakistan. Saboori et al. [65] for Malaysia and Shahbaz et al. [68] 
for Romania. 

The short run dynamics are also reported in lower segment of 
Table 6 and the results indicated that coal consumption is also 
main contributor to CO, emissions and it is statistically significant 
at 1% level. Interestingly, international trade reduces energy 
emissions significantly at 5% significance level. The coefficients 
of linear term of real GDP per capita and nonlinear i.e., squared 
term of real GDP per capita are positive and negative, respec- 
tively. These estimates are significant at 1% significance level. This 
further confirms the existence of EKC in case of India. The short 
run estimates are less than long run estimates which indicates the 
reliability and stability of estimated results. 

The estimated coefficient of the lagged ECM term is — 0.2994 and 
significant at the 1% level. This establishes long run relationship 
among the running variables. This suggests that changes in CO2 
emissions from the short run to the long run are corrected by 29.94% 
each year. The significance of lagged error correction term further 
confirms our established long run relationship between coal con- 
sumption, economic growth, trade openness and CO, emissions. 

The diagnostics tests such as LM test for serial correlation, 
normality of residual term and white heteroscedisticity test for 
short run model are reported in Table 5. The findings suggest that 


Fig. 3. Plot of cumulative sum of recursive residuals. The straight lines 
represent critical bounds at 5% significance level. 


the short-run model passes all diagnostic tests successfully. The 
evidence indicates that error term is normally distributed and 
there is an absence of serial correlation. There is no evidence of 
autoregressive conditional heteroscedisticity and same inference 
for white heteroscedisticity. The short run model has passed 
Remsay reset test which confirms that the functional form of 
short run model is correctly specified. 

It is suggested by Browns et al. [53] that stability of long run 
and short run parameters can be tested by applying CUSUM and 
CUSUMszg tests. In doing so, the cumulative sum (CUSUM) and 
cumulative sum of squares (CUSUMsq) tests have been employed 
for parameter stability. 

The results are shown in Figs. 3 and 4. The graph should lie 
within the appropriate bounds if the parameters are 
stable (Pesaran et al. [51]). The graph of CUSUM at the 5% level 
and we can claim the stability of long-and-short runs parameters. 


4.1. The VECM Granger causality analysis 


The presence of cointegration among the variables implies that 
causality relation must be existed at least from one side. The 
directional relationship between coal consumption, economic 
growth, international trade and CO, emissions will provide help 
in articulating comprehensive policy to economic growth by 
controlling environment from degradation and utilize energy 
efficient technologies importing from advanced countries. We 
applied Granger causality test within the VECM framework to 
detect the causality between the variables. Table 6 reports the 
results of the VECM Granger causality analysis. The long run 
causality is captured by a significant t-test on a negative coeffi- 
cient of the lagged error-correction term ECM;_;. The jointly 
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1985 1990 1995 2000 2005 2010 


Fig. 4. Plot of cumulative sum of squares of recursive residuals. The straight 
lines represent critical bounds at 5% significance level. 


significant LR test on the lagged explanatory variables shows 
short-run causality. 

In long run, our empirical evidence reports bidirectional 
causality relationship between income and CO2 emissions. The 
feedback hypothesis exists between coal consumption and CO, 
emissions. Coal consumption and economic growth Granger cause 
each other. This finding is contradictory with Li and Li [62] who 
reported that coal consumption Granger causes economic growth 
while consistent with Paul and Bhattacharya [63] showed feed- 
back between energy consumption and economic growth. This 
suggests that government should encourage energy exploration 
polices which in resulting spur economic growth and hence 
demand for energy. There is bidirectional causal relationship 
between economic growth and CO2 emissions. These findings 
confirm our long run as well short run results that income has 
positive and negative effects on CO emissions. This further 
validates the existence of environmental Kuznets curve (EKC). 
On contrary, Ghosh [26] and Alam et al. [60] presented no causal 
relationship between income and CO: emissions. 

This shows the significance of our empirical evidence as we 
applied more specified economic model. Trade openness Granger 
causes economic growth, coal consumption and CO, emissions. 
This implies that trade openness stimulates economic activity and 
adds in energy demand (coal consumption) which emits CO2 
emissions. In short run, bidirectional causality is found between 
economic growth and CO2 emissions, coal consumption and CO2 
emissions, trade openness and CO, emissions. The findings 
indicated that the long run income elasticity for CO2 emissions 
is less than the short run elasticity is further evidence for an EKC 
(see Narayan and Narayan [30] for more). The feedback hypoth- 
esis also exists between trade and economic growth. The unidir- 
ectional causal relationship is found running from trade openness 
to coal consumption and no causality exists between coal con- 
sumption and economic growth. The results of joint causality 
confirm the long-and-short runs findings. 


5. Conclusion and policy implications 


This paper investigated the dynamic linkages between coal 
consumption, economic growth, trade openness and CO, emis- 
sions in case of India over the period of 1966-2011. We have 
applied Ng-Perron [57] unit root test to test stationarity of the 
variables and Narayan and Popp [49] is accompanied to examine 
the robustness of empirical evidence. Further, the ARDL bounds 
testing approach to cointegration has employed to investigated 
the long run relationship between the variables and Johansen 
multivariate cointegration is applied to test the strength of long 
run results. The present study has important policy implications 
for Indian economy. The empirical evidence of strong relationship 
between coal consumption and economic growth indicates that 


sustained level of economic growth can be achieved by exploring 
the new sources of energy supply for long span of time. 

Our results confirmed that all the variables are cointegrated 
for long run relationship and long run results are robust. Coal 
consumption is main contributor to CO, emissions. A rise in 
economic growth raises demand for coal that emits energy 
pollutants. Trade openness is also positively linked with CO, 
emissions but in short run former declines later. Results also 
validated the existence of environmental Kuznets curve (EKC) in 
long run as well as in short run in case of India. The causality 
analysis explored the bidirectional causal relationship between 
economic growth and CO, emissions and feedback hypothesis 
exists between coal consumption and CO% emissions. Coal con- 
sumption and economic growth Granger cause each other. Trade 
openness Granger causes economic growth, coal consumption 
and CO2 emissions. Economic growth and CO2 emissions Granger 
cause trade openness in short run. 

Hence, our results shows that dependence on the coal con- 
sumption as a factor that promotes economic growth might be 
feasible option for short run but in the long run we have to look 
for alternative of this. This is particularly because coal consump- 
tion increases CO, emissions and which indirectly has negative 
impact on the economic growth. The other interesting issues are 
the Granger-causal relationship between CO, emissions, trade 
openness, and coal consumption. Since trade openness is found to 
Granger-cause CO z emissions, we can suggest Indian policy- 
makers on a few important points: (1) allowance of foreign 
investors to invest in the India particularly in know-how and 
advanced technologies (2) free entry of FDI because when FDI 
flows into a host country, FDI may transfer new ideas, technol- 
ogies and managerial skills to domestic firms and which in turn 
expected to reduce CO, emissions (3) promotion of such activities 
through which foreign firms may provide continuous technologi- 
cal assistance, such as introduction of new technologies and 
training of workers, to the domestic firms which will lead to 
decrease in CO, emission in one hand and increase economic 
growth on the other. These processes are also expected to reduce 
dependence on the consumption of coal without harming the 
interest of higher economic growth. 
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